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EXPERIMENTALLY DETERMINED STIFFNESS AND DAMPING OF AN 
INHERENTLY COMPENSATED AIR SQUEEZE-FILM DAMPER 
by Robert E. Cunningham 
Lewis Research Center 

SUMMARY 

Experiments were conducted to determine the dynamic characteristics of damping 
and stiffness of an externally pressurized, inherently compensated, compressible 
squeeze-film damper. Each damper was 76. 2 millimeters (3 in. ) square and contained 
20 equally spaced feed holes 1. 5 millimeters (0. 059 in. ) in diameter. Results were ob- 
tained at clearances of 63.5, 50.8, 38.1, and 25 micrometers (0.0025, 0.002, 0.0015, 
and 0.001 in. ), at supply pressure ratios of 2 and 3, and for excitation frequencies to 
36 000 cycles per minute. The test apparatus consisted of three such dampers arranged 
in an equilateral, triangular pattern. 

Experimental values of the squeeze-film damping when compared with a small pres- 
sure perturbation analysis were much greater than predicted at low squeeze numbers 
(<20) and slightly less than predicted at higher squeeze numbers (>20). 

It was determined for limited applications requiring only light damping in a flexible 
bearing support that a gas squeeze -film damper similar to the type investigated could 
supply approximately the same amount of damping as cotild the previously designed oil 
squeeze -film damper. The clearances in the gas squeeze -film damper would have to be 
much smaller to obtain the same magnitude of damping - that is, 38. 1 micrometers 
(0.0015 in. ) for the air as compared to 127 micrometers (0.005 in. ) for the oil. 

Experimental stiffness values when compared to plotted results of the small pertur- 
bation analysis were much less than predicted for low squeeze numbers and much greater 
for high squeeze numbers. 


INTRODUCTION 

Small, lightweight, high speed turbomachines are being designed in increasing 
numbers. To insure the safe operation of these machines and to provide reasonable 



operating longevity, considerable attention must be given to the interaction of the rotor 
and its bearing supports. 

The rotors of these machines can be designed to operate either as rigid or flexible 
members. The bearing supports can be an integral part of the machine frame, having a 
stiffness of the same order of magnitude as the rolling element bearings, or they can be 
flexible, having a stiffness much less than that of the bearings. These so-called ’’soft” 
supports serve two functions: 

(1) They can suppress the system rigid body critical speeds well below the operating 
range. 

(2) They can attenuate the forces transmitted to the bearii^s (ref. 1). 

Currently , a number of production turbojet engines are using a type of cantilevered 
’’squirrel cage” spring to support the bearings. Once support flexibility is introduced, 
however, it becomes necessary to also provide a damping mechanism. Flexibility with- 
out damping can produce a system highly sensitive to shock as in the case of suddenly 
applied xmbalance loads (ref. 1). 

The literature on rotor dynamics contains many reports which treat the dynamic 
response of a rigid rotor in rigid supports. References 2 and 3 treat the more complex 
problem of an elastic rotor in elastic, damped supports. With the aid of the computer 
programs developed in these analyses, the machinery designer can, with some degree of 
confidence, predict the response of an unbalanced rotor to assigned values of support 
stiffness and damping. To build into the machine just the correct amoimt of stiffness 
and damping is, however, a much more difficult task. 

More often than not, compromise solutions are necessary to satisfy the various sets 
of operating conditions encountered in high speed machinery. For example, to guard 
against bearing failures due to large transient overloads, much higher damping in the 
supports is used than is necessary to operate through a critical speed at nominal unbal- 
ance conditions. The greater damping can, at operating speeds, transmit larger forces 
to the bearings and thus result in their premature failure (ref. 1). 

Some of the more conventional methods of obtaining support damping (lubricating 
oils and elastomers) have reduced effectiveness when exposed to high temperatures in an 
actual machine. Elastomers used as springs and/or dampers can have their d 3 mamic 
properties significantly altered by high temperatures, either self -induced or from adja- 
cent hot parts. 

The high temperatures in bearii^ compartments can change the viscosity of lubri- 
cating oil and thus result in a change in the damping coefficient. In certain applications 
the temperatures may be high enough to even prohibit the use of such oils as damper 
fluids. 

For these reasons one is led to consider an externally pressurized, air squeeze - 
film bearing to provide the flexibility and damping needed to attenuate the forces trans- 
mitted to the supports in certain high temperature applications. 
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The gas bearing literature contains a number of reports which investigate the char- 
acteristics of hydrodynamic, compressible, squeeze-film bearings (refs. 4 and 5). 

There are also several reports on the dynamic characteristics of the externally pres- 
surized bearii^ (refs. 6 to 8). Reference 9 treats the dynamic characteristics of an in- 
herently compensated, multiple feed hole, circular gas thrust bearing. The research 
reported in references 7 and 10 includes experimental data on the dynamic properties of 
externally pressurized joiornal bearings. 

To be able to use such a bearing as a flexible support, the stiffness and damping 
properties must be predictable. To date, only a small amount of experimental data is 
available to compare with the analyses. One objective of this investigation is to deter- 
mine the stiffness and damping of an externally pressurized, inherently compensated, 
air film damper at clearances of 63.5, 50. 8, 38. 1, and 25 micrometers (0. 0025, 0. 002, 
0. 0015, and 0. 001 in. ) for supply pressure ratios of 2 and 3 and excitation frequencies 
to 36 000 cycles per minute. The results are compared to a small pressure perturbation 
linearized analysis for a multiple feed hole, inherently compensated, gas film damper 
(ref. 11). 

A second objective of this investigation is to determine the feasibility of usii^ an 
arrangement of three externally pressurized bearings as a flexible damped support 
whose dynamic characteristics can be readily altered to attenuate amplitudes of motion 
and transmitted forces while the rotor operates through its critical speeds. 

The polygonal arrangement of the three rectangular shaped, externally pressurized 
pads used as a flexible, damped support lends itself to ease of clearance adjustment. 

This would be very difficult to do with a circular bearing. The rectangular pads ar- 
ranged in a triangle also prevent rotation without using anti -rotation pins, etc. Finally, 
the pressure in each individual pad can be regulated to obtain the amo\mt of stiffness 
necessary to satisfy the rotor unbalance forces at any particular speed. This would also 
be difficvilt to do in a support with circular geometry. Because of the unique features 
embodied in this concept, a patent disclosure has been filed. 


SYMBOL LIST 

A ratio of center region length to total pad length, Z/L 
B squeeze film dampir^ coefficient, N-sec/m ((lb-sec )/in. ) 
B dimensionless damping, Bh /juL 
Cj^ orifice discharge coefficient 
d^ orifice diameter, mm (in.) 
e housing peak to peak amplitude, mm (in. ) 

F rotating unbalance force , N (lb) 
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f coefficient representii^ integral of pressure gradient along source boundary, 

dimensionless 

9 p 

g acceleration of gravity, m/sec (in. /sec ) 
h film thickness, mm (in. ) 

K squeeze film stiffness, N/m, (Ib/in.) 

K dimensionless stiffness, Kh/(L^Pg^[(Pg/Pg^) - 1 } 
k adiabatic gas exponent 

I center region length, m (in. ) 

L damper pad length, m (in.) 

M housing mass, kg (lb-sec /in.) 
m unbalance mass, Isg (lb-sec /in. ) 

N turbine speed, rev /min 

n number of orifices 

p p 

P^ ambient pressure, N/m (Ib/in, abs) 

2 2 

P„ supply pressure, N/m (Ib/in. abs) 

R gas constant, J/(kg)(K) (in. -lb/(lb)(°R)) 

r distance to unbalance weight, m (in. ) 

T temperature, K 

t time, sec 

U unbalance, kg-m (oz-in. ) 

W housing weight, kg (lb) 

w unbalance weight, kg (lb) 

X, y, z coordinate axes 

A restrictor coefficient, 6p C^nyrd^ |2kRT/(k - 1) 

p absolute viscosity, N-sec/m^ (Ib-sec/in. ^) 

2 2 

(7 squeeze number, 12pojL /h P 

Ci 

(p phase angle, deg 

o) angular velocity, rad/sec 

Wp angular velocity at resonance, rad/sec 

Subscript: 
e error value 


l/2]/p^h2f 
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APPARATUS 


Mechanical Components 

The experimental results presented in this report were obtained using the test appa- 
ratus shown schematically in figxire 1 and photographically in figures 2 and 3. An equi- 
lateral triangular damper housing having three precision ground and lapped faces is lo- 
cated by three externally pressurized gas damper pads. These pads are 7. 62 milli- 
meters (3 in. ) square and contain twenty 1. 50 -millimeter- (0. 059 -in. -) diameter feed 
holes arranged in a square pattern (fig. 4). Two of the pads are rigidly attached to a 
base plate while the third is mounted on a micrometer slide permitting adjustment of the 
gas film clearance in all three damper pads. Each pad has its own air supply which can 
be regulated independently. This pressure was measured by Bourdon tube type gages. 

The housing was supported vertically by a pocketed externally pressurized thrust bearing. 

An air driven tm-bine wheel is moimted on the inside of a circular cavity in the hous- 
ing whose plane is perpendicular to the pad faces (see fig. 1). The turbine rotates on a 
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ball bearing and contains a small unbalance moment of 3.99x10 kilogram -meter 
(0. 0553 oz-in. ). This amount of unbalance was measured on a two plane dynamic balanc- 
ing machine. The rotating timbine wheel produces a harmonic exciting force acting nor- 
mal to each gas bearing surface; this is shown schematically in figure 5. 


INSTRUMENTATION 

The instrumentation used to obtain the experimental results is shown schematically 
in figvme 6. Two noncontacting, proximity measuring probes are mounted flush with the 
surfaces of each of the three pads. The output voltage from the probes is proportional to 
this displacement. This signal is amplified, fed into a tracking filter from where it goes 
to a peak to peak measuring voltmeter, and finally displayed on an X-Y-Y' plotter. 

A magnetic probe mounted over the barbine wheel produces a conditioned voltage 
pulse once per revolution when triggered by a single narrow tooth on the turbine wheel 
located at the same angle as the unbalance weight. This pulse, in addition to measuring 
the angular velocity of the turbine wheel, provides the reference for measiiring phase 
angle. This is the angle between the force and displacement vectors as measured by an 
electronic phase meter. The output from the phase meter is then fed to the Y*-axts of 
the X-Y-Y' plotter. A dynamic input sine converter (DISC) takes the repetitive pulse and 
converts it to a dc voltage which is proportional to frequency. This signal is then ap- 
plied to the X-axis of the plotter where housing amplitude and phase angle are automati- 
cally plotted against rotor frequency. 
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The DISC is also used to convert the dc pulse from the magnetic pickup to a sinus- 
oidally varying output which in tiirn is used to drive the carrier generator of the tracking 
filter. The tracking filter has both a 10- and 20 -hertz bandpass filter, either one of 
which can be used. Amplitude and phase can also be measured on a cathode ray oscil- 
loscope (CRO). The reference pulse relative to the sinusoidal motion of the housing is 
shown in figure 7 ; the photographs of the CRO screen were taken during a test run. 


PROCEDURE 
Initial Assembly 

In order to achieve the best alinement possible between the damper faces and the 
triangular housing, the following steps were taken in the fabrication and initial assembly 
of these components. 

The three faces of the damper housing were ground and lapped to within ±5 micro- 
meters (±0. 0002 in. ) of a true equilateral triangle. In addition, each face was held with- 
in three helium light bands of being optically flat. These same tolerances were also held 
on a circular externally pressurized thrust bearing and on the mating surface at the bot- 
tom of the triangular damper housing. 

Before attaching any of three damper pads to the base plate, the pads were first 
moved into contact with the triangular housing after having positioned 10-micrometer- 
(0. 0004-in, -) thick paper shims at each of the four corners of each pad. Shims were 
placed, if required, between the base plate and the base of each externally pressurized 
pad until equal tension was achieved on each of the paper shims. The pads were then 
carefully bolted down and doweled in place. The same procedure was used to insure the 
parallelism of the adjustable pad with its matir^ svmface on the housing when bolted to 
the micrometer adjustable slide. 


PRETEST 

Before each test the clearance between each bearing and its respective mating sur- 
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face was adjusted. First, air was supplied to the thrust bearing at 55. 2x10 newtons per 
square meter (80 psig in the plenum) to lift the housing from the thrust bearing surface. 
When the housing floated freely at approximately a 38. 1-micrometer (0.0015-in.) clear- 
ance, the adjustable pad was moved in until metal to metal contact was made on all three 
pads. The two capacitance probes in each pad were adjusted to approximately 50. 8 mi- 
crometers (0. 002 in. ) below the surface of each pad. This was determined by the output 
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voltage from the probes as read on a digital voltmeter. This value was recorded and be- 
came the zero clearance reading. 

The adjustable pad was then moved back with the micrometer slide until a specific 
clearance was achieved. The air supply was then turned on to each pad and the pressime 
in each plenum adjusted to the same value. The clearance in each pad was checked on 
the voltmeter and recorded. Stops on the adjustable damper pad were then positioned so 
that the play in the micrometer slide did not allow any movement of this pad. The aver- 
age clearance was within 5 percent for all three pads. 

The next steps involved calibrating the instrumentation shown in figure 6. Prior to 
each test the two tracking channel bandpass filters of 10 and 20 hertz were calibrated us- 
ing an audio oscillator monitored by an electronic counter. The amplitude and phase 
axes were then calibrated using the same audio oscillator, dynamic sine converter, and 
dc peak to peak reading voltmeter. 


TEST 

The air supply valve on a speed controller was opened to provide drive air to the un- 
balanced turbine wheel. The plotter was turned on to record the amplitude and phase of 
the resulting housing motion relative to a single pad over a range of frequencies to 
36 000 cycles per minute. A cathode ray oscilloscope was also used to monitor both the 
amplitude and phase as shown in figure 7. The test clearance and zero clearance read- 
ings were again checked at the conclusion of each test. 


RESULTS AND DISCUSSION 

The experimental results of this investigation are presented in figures 8 to 15 and 
are compared to the theoretical results obtained in the small pressure perturbation lin- 
earized analysis of reference 11. 

Experimental values of damping and stiffness were obtained from the equation of 
motion for a damped single degree of freedom spring mass system subjected to a har- 
monic exciting force. The following force displacement equation is from reference 12 
and pertains to the system under investigation: 

M jeo)^ sin (wt - - — j^ew sin ^cot - <p + — ^ j^e sin (cot - + F sin cot = 0 (1 ) 


From this equation these vector diagrams can be constructed. 






Damping and stiffness coefficients are solved for as follows using the vector diagram on 
the right: 

Damping coefficient: 


B = 


2F sin <p 
3eco 


N-sec 

M 



Stiffness coefficient: 


_ 2F cos (p 
3e 


^ 2Mct)^ 

3 



( 2 ) 


(3) 


The measured variables were peak to peak amplitude, phase angle <p, and frequency co 
of the excitir^ force F which was equal to the turbine frequency. 

Since there are three pads, there are actually three equivalent springs and dampers. 
However, it is easily verified that when all pads are identical the spring and dampii^ 
constants in any direction are equal to (3/2)Kand (3/2)B, respectively, where K and 
B are constants for a single pad. Since the constants are directionally uniform, a ro- 
tating load, as applied in the apparatus, produces a circular motion. This was verified 
by accelerometers mounted orthogonally on the damper housii^. The accelerometer sig- 
nal, displayed on an X-Y curve tracing oscilloscope, showed a circular orbit. 

Figure 8 is for an actual plot of the housing peak to peak motion and the phase angle 
as a fimction of excitation frequency. Figure 8(a) is at a clearance of 63.5 micrometers 
(0. 0025 in. ) and shows a large buildup in amplitude at the resonant frequency; this indi- 
cates very light damping. Figure 8(b) is at a clearance of 25 micrometers (0. 001 in.) 
where no peak amplitude occurs at resonance; this indicates the presence of much 
greater damping. 

In order to compare the experimental values of damping and stiffness with the theo- 
retical values as plotted in reference 11, the data were put in nondimensional form by the 
following equations: 

Dimensionless damping: 
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(4) 


B = 


Bh' 

fiL 


Dimensionless stiffness; 


Kh 


(5) 


Squeeze number: 


a 


_ 12]ua)L' 
h^P„ 


( 6 ) 


DAMPING 

Figxire 9 is for a pressure ratio of 2. The theoretical dimensionless dampii^ re- 
mains relatively constant through squeeze numbers of approximately 20; beyond 20 there 
is a pronounced reduction in dampii^. For very high squeeze numbers the predicted 
damping diminishes to extremely small values (ref. 11). At an initial clearance of 63. 5 
micrometers (0. 0025 in.) the experimental data from six tests (fig. 9(a)) when compared 
with theory show that the damping is much higher than predicted. In addition, the data 
show considerable scatter. This is probably due to the fact that the initial clearance h 
was reset between tests and that tests at the same clearance were not run consecutively. 
For successively smaller initial clearances, figures 9(b) to (d) show fair agreement with 
theory. At values of 38. 1 and 25 micrometers (0. 0015 and 0. 001 in. ) (figs. 9(c) and (d)), 
the theory overpredicts throughout the range of squeeze numbers Investigated. Similarly, 
at a supply pressure ratio of 3 (figs. 10(c) to (d)) the experimental results exhibit the 
same trend with slightly better agreement for clearances of 38. 1 and 25 micrometers 
(0.0015 and 0.001 in.) (figs. 10(c) and (d)). 

Figures 11 and 12 are cross plots of the dimensional damping coefficient B against 
initial film clearance at a squeeze number ct of 10 and supply pressure ratios P_/Pa 
of 2 and 3, respectively. A squeeze number of 10 was chosen for this comparison be- 
cause the damping data were taken at this value for all values of clearances tested. 

These two figures show how rapidly the squeeze film damping increases with decreasii^ 
clearance. The actual measured damping is less than predicted at clearances up to 
50.8 micrometers (0.0020 in.) at P„/P„ = 2 and up to 45.7 micrometers (0.0018 in.) 
at P /P = 3; the actual damping becomes greater than predicted at clearances larger 
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than these. These damping curves are typical of externally pressurized gas bearings, 
and similar results are shown in reference 7. 

STIFFNESS 

Gas film stiffness in nondimensional form is plotted as a function of squeeze number 
in figures 13 and 14. The experimental values are compared with the theory of refer- 
ence 11. Agreement between the experimental data and the theory is poor at the higher 
clearances. For a clearance of 25 micrometers (0. 001 in.) the agreement is fair at sup- 
ply pressure ratios of 2 and 3. Values of stiffness increase at a much more rapid rate 
for increasing squeeze numbers than the theory predicts. At low squeeze numbers the 
measured stiffness is much less than theory predicts except at a clearance of 25 micro- 
meters (0. 001 in. ) and a supply pressure ratio of 2. 

In figure 15 theoretical and measured stiffnesses (figs. 13 and 14) are plotted as 

functions of the initial clearance at a squeeze number of 10 and for supply pressure ratios 

of 2 and 3. Experimental values of squeeze film stiffnesses ranged from a high of 
5 

143.7x10 newtons per meter (82 000 Ib/in. ) at a supply pressure ratio of 3 and a clear- 
ance of 63. 5 micrometers (0. 0025 in. ) to a low of 13. ixlO newtons per meter (7500 lb/ 
in. ) at a supply pressure ratio of 2 and a clearance of 25 micrometers (0. 001 in. ). 


ERROR ANALYSIS 

The dynamic characteristics of a gas film, specifically stiffness and dampii^, can- 
not be measured directly. Values of these coefficients must be calculated from the 
values that can be measured directly - these being displacement, phase angle, frequency, 
and initial clearance. 

A number of different sensors and instruments are necessary to measure these var- 
iables and, of course, each is subject to error. The following are typical accuracies of 
instruments used as obtained from manufacturers' operating maniials: 



Percent error 

Proximity sensor output 

±2.6 

Digital voltmeter 

±3.0 

Amplifier 

±.01 

Dynamic input sine converter 

CM 

-H 

Tracking filter 

±2.0 

Phase meter 

.08 

Peak to peak dc voltmeter 

o 

±1 

X-Y-Y‘ plotter 

±.2 
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In addition to the errors inherent in the instrumentation, there is a probable misaline- 
ment error of ±0. 0002 inch. 

The following sample calcvilations illustrate the potential error in the measured 
values of amplitude of housing motion, frequency of exciting force, the phase angle, and 
initial clearance setting. The measured values were at test conditions of P_/P„ = 2 

O <4# 

and h = 63. 5 micrometers (0. 0025 in. ). 

(1) Potential error in peak to peak amplitude - dependent on proximity sensor, am- 
plifier, tracking filter, peak to peak voltmeter, and X’-Y-Y' plotter (observed value, 

e = 13.46 pm (0.00053 in.)): 

eg = 13. 46 pm (0. 00053 in. ) X ±0. 088 = ±1. 18 pm (4.63x10"^ in. ) 
eg = 13.36 ± 1.18 = 14. 5 to 12.2 pm (0. 572x10"^ to 0.48x10"^ in.) 

(2) Potential error in frequency (signal from sine converter to X-Y-Y’ plotter (ob- 
served value, N = 24 000 cycles/min): 

Ng = 24 000 cycles/min x ±0. 0045 = ±108 cycles/min 

Ng = 24 108 to 23 892 cycles/min 

(3) Potential error in phase angle, from phase meter to X-Y-Y' plotter (observed 
phase angle, 171°): 


<Pg = 171° X ±0.0028 = ±0.479° 

(^g = 171.5° to 170.5° 

(4) Potential error in initial clearance setting from proximity sensor and digital 
voltmeter (observed clearance, h = 63.5 pm (0.0025 in.)); 

hg = 63. 5 pm X ±0. 0056 = ±0. 356 pm (±1. 4X10"® in. ) 

hg = 63. 86 to 63. 14 pm (0. 00251 to 0. 00249 in. ) 


To this value must be added the possible misalinement error of ±5. 08 pm (±0. 0002 in. ): 
hg = 63.86 ± 5.08 = 68.94 to 58.78 pm (0.00271 to 0.00231 in.) 
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Taking a combination of measured values that would give the highest and lowest pos- 
sible errors results in dimensionless damping values that could be 49 percent too high or 
31 percent too low. 

The following sample calculations illustrate the variations in the experimental values 
of damping, stiffness, and squeeze numbers; 

Damping coefficient, B; 


g ^ 2. 79x10 N sin <p 
e 

R = Rvin-8 Nsin cp 

e 

eg = 14. 5 to 12.2 pm (0.572x10'^ to 0.48x10“^ in.) 

Ng = 24 108 to 23 892 cycles/min 
<Pg = 171.5° to 170.5° 

hg = 68. 9 to 58.7 pm (0.00271 to 0.00231 in.) 

Substituting the positive and negative error values into the equation for damping coeffi- 
cient gives 


N-sec 

m 

lb -sec 
in. 


Bg = 908 to 676 N-sec/m (5. 18 to 3. 86 Ib-sec/in. ) 

as compared to the nominal value of 776 newton-seconds per meter (4.43 lb -sec/in.). 
Dimensionless damping, B: 


B = 


Bh 


3 


pL 


4 


p = 1.81x10"^ N-sec/m^ (2. 62x10'^ Ib-sec/in.^) 
L = 7.62 mm (3 in.) 

Substituting the values into the equation gives 


Bg = 0.486 to 0.224 

as compared to the nominal value of 0,326. 

Squeeze number, a; 
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^ _ 12jua;L^ _ 2.96 x10~^N 

Ng = 24 108 to 23 892 cycles/min 

hg = 68.9 to 58.7 jim (0.00271 to 0.00231 in.) 

P, = 98 289 N/m^ (14. 25 psia) 

d. 

Substituting the values gives 

CT =9.38 to 6.76 

as compared to the nominal value of 7. 98. 

Instrumentation errors in the experimental values of stiffness could result in values 
that are greater by 22 percent and lower by 23. 8 percent. Error values of measured am- 
plitude, phase angle, frequency, and initial clearance are given again as follows: 

eg = 19. 5 to 12. 2 jum (0. 572x10“^ to 0. 48x10"^ in. ) 

Ng = 24 108 to 23 892 cycles/min 

cp^ = 171. 5° to 170.5° 

hg = 68. 9 to 58. 7 jum (0. 00271 to 0. 00231 in. ) 

Substituting these values into the expression for stiffness coefficient, 

^ _ 2(Uu>h cos <p ^ 2Mo? N 

3e 3 m 

results in values of 

Kg = 87. 3X10® to 63. 7x10® N/m (49 843 to 36 362 Ib/in. ) 

as compared to the nominal value of 77. 3x10® newtons per meter (44 142 Ib/in. ). 
Dimensionless stiffness is defined as 



13 


Kh 


K = 

l] 

Substituting the values of error values of into this equation gives values of dimen- 
sionless stiffness of 


Kg = 1.05 to 0.655 

as compared to the nominal value of 0. 86. The following, taken from reference 13, is 
offered as an explanation for the disagreement between the experimental values of damp- 
ing and those predicted by the theory of reference 11. For a self-acting gas lubricated 
squeeze -film bearing, the film is effectively incompressible when the normal velocities 
are relatively low, that is, for low squeeze numbers. This means that the viscous or 
damping forces predominate. At high relative velocities the film approaches a purely 
compressible state where there is virtually no fluid flow. This is, however, an exter- 
nally pressurized bearing. There is always gas flow out over the sill (regardless of 
squeeze number), and, therefore, the measured damping does not decrease as rapidly as 
predicted. The analysis of reference 11 is based on a small perturbation linearized so- 
lution. This assumption, however, would not account for the disagreement since the 
maximum eccentricity ratio observed in these tests was 0. 3. 

The explanation given earlier in this discussion for the variance observed between 
experimental values of damping applies to the film stiffness also (ref. 13): When the 
velocity of the moving surface is relatively low and clearances are high (low squeeze 
number, a — 0), the film is effectively incompressible with high viscous forces. This, 
of course, applies for a pure h 5 ^drodynamic bearing and is not completely true for an 
externally pressurized bearing of the type being considered. 

At very small clearances and high approach velocities a somewhat compressible 
film state exists where the hydrodynamic component of the flow is reduced. This con- 
dition results in much greater stiffness values than the theory is capable of predicting. 

For the smaller clearance and a supply pressure ratio of 2 there is good agreement 
between theory and experiment. The measured stiffness, however, continues to in- 
crease over the range of clearances investigated (25 to 63.5 pm (0. 001 to 0. 0025 in. )) 
whereas the predicted stiffness peaks and then drops off. Where the theory shows a sig- 
nificant increase in film stiffness when the supply pressure ratio is increased from 2 
to 3, the measured stiffness increase is only moderate. 
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COMPARISON BETWEEN OIL AND GAS FILM DAMPER 


One of the objectives of this study was to determine the feasibility of using an inher- 
ently compensated, externally pressurized, gas film bearing as a flexible damped sup- 
port. The test apparatus was designed to provide the capability of readily changing both 
the stiffness and damping characteristics without shutting down the drive system. Con- 
trol was provided of both the supply air to the pads and the initial clearance by a micro- 
meter adjustment. 

Figures 11 and 12 show the variations, both experimental and predicted, in damping 
coefficients of a single pad for a range of clearances. Average values of the measured 
damping for a single pad range from a high of 3294 newton -seconds per meter (18. 8 Ib- 
sec/in. ) at a supply pressure ratio of 3 and a clearance of 25 micrometers (0. 001 in. ) 
down to 561 newton -seconds per meter (3, 2 lb -sec/in. ) at a clearance of 25 micrometers 
(0. 0025 in. ) and a supply pressure ratio of 2. 

It was suggested in the INTRODUCTION that an externally pressurized gas film 
damper of the type investigated might find application in small, lightweight turbocom- 
pressors. It could be used where bearing compartment temperatures were so high so as 
to preclude using either oil or elastomers in energy dissipation devices. The calcula- 
tions that follow show that the type of gas film damper tested can provide sufficient damp- 
ing to malce it a candidate for replacii^ the oil squeeze -film damper. The oil damper 
dimensions and the amount of support damping required were calculated for a small, 
lightweight turborotor simulator (described in ref. 14). This simulator is required to 

_4 

operate above its first bending critical speed with an imbalance of 2. 88x10 kilogram - 
meter (0. 4 oz-in. ). 

The rotor weighs 5. 68 kilograms (12. 52 lb) and is supported in two ball bearings. 

This amount of unbalance is large for a rotor of this weight; however, it represents an 
unbalance condition that could result if the original balance changed because of blade 
erosion or wear of machine components, 

The equation for dimensionless damping of an Incompressible lubricant was obtained 
from reference 15: 


B = Bc|/2)ul| 

where 

B dimensionless damping for oil damper at e = 0.4 
L damper length, 0.01 m (0.45 in.) 

R damper radius, 0.04 m (1.56 in.) 
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ju absolute oil viscosity, 0. 12 N-sec/m^ (1. 73x10 Ib-sec/in, 

Cj^ damper clearance, 127 pm (O.OOSin. ) 

Substituting these values into the previous expression and solving for the damping coeffi- 
cient give a B value of 1405 newton-seconds per meter (8.02 lb -sec/in. ). 

A supply pressure ratio P„/P„ of 3, a squeeze number of 10, and a clearance of 

S d- 

38 micrometers (0. 0015 in. ) give a damping coefficient of 1174 newton-seconds per meter 
(6. 7 Ib-sec/in. ) as shown in figure 12. This value of the damping coefficient is for a sin- 
gle pad only. If three pads are used as arranged in figure 5, the total damping coeffi- 
cient would be 1. 5x1174 or 1840 newton-seconds per meter (10. 05 Ib-sec/in.). Thus, 
the air film damper would provide sufficient damping for the application although as 
presently designed the damper is larger and the clearances are much smaller. 


SUMMARY OF RESULTS 

Experimental values of damping and stifftiess for an externally pressurized, inher- 
ently compensated, air film damper were determined for clearances of 63. 5, 50. 8, 38. 1, 
and 25 micrometers (0. 0025, 0. 002, 0. 0015, and 0. 001 in. ) at frequencies up to 36 000 
cycles per minute and supply pressure ratios of 2 and 3. The damper was 76.2 mil- 
limeters (3 in.) square and contained 20 equally spaced feed holes 1.4 millimeters 
(0. 059 in. ) in diameter. The test apparatus had three such pads arrai^ed in an eqm- 
lateral triangular pattern so that the air film in each pad could be excited by a rotating 
unbalance force. 

Measimed squeeze film damping coefficients for a single bearii^ ranged from a high 
of 3294 newton-seconds per meter (18.8 Ib-sec/in.) at a clearance of 25 micrometers 
(0. 001 in. ) and a supply pressure ratio of 3 down to 67. 1 newton -seconds per meter 
(3. 2 Ib-sec/in. ) at a clearance of 63. 5 micrometers (0. 0025 in. ) and a supply pressure 
ratio of 2. 

The experimental damping values at a supply pressure ratio of 2, clearances of 
63.5 and 50.8 micrometers (0.0025 and 0.002 in.), and squeeze numbers less than 20 are 
considerably higher than those predicted by a small pressure perturbation theory. At 
clearances of 38.1 and 25 micrometers (0.0015 and 0.001 in.) the measured damping is 
less than that predicted by theory. At a supply pressure ratio of 3 the measured damp- 
ing is again greater than that predicted at high clearances and squeeze numbers less 
than 20. Agreement is significantly better, however, at clearances of 38.1 and 25 mi- 
crometers (0.0015 and 0.001 in.) and relatively high squeeze numbers. 

Measured values of the squeeze-film stiffness ranged from a low of 13. ixlO 
newtons per meter (7500 Ib/in. ) at a clearance of 25 micrometers (0. 001 in. ) and a sup- 
ply pressure ratio of 2 to a high of 143.7x10^ newtons per meter (82 000 Ib/in. ) at a 
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clearance of 63. 5 micrometers (0. 0025 in. ) and at a supply pressure ratio of 3. 

Measured values of the stiffness coefficient are generally in poor agreement with the 
theoretical values, particularly at high clearances and low squeeze numbers; that is, 
they are much less than the predicted values. As was the case with the damping coeffi- 
cients, the measured stiffnesses at small clearances and high squeeze numbers showed 
significantly better agreement although they were slightly higher than predicted. 

The present damper design provides a damping coefficient of 1840 newton-seconds 
per meter (10. 05 Ib-sec/in. ) for a pressure ratio of 3, a clearance of 38 micrometers 
(0.0015 in. ), and a squeeze number of 10. This compares with 1405 newton -seconds per 
meter (8. 02 Ib-sec/in. ) for a previously designed oil squeeze -film damper with a sim- 
ilar diameter. 

Results of an error analysis showed that values of the experimental dimensionless 
damping could be 49 percent higher or 31 percent lower than the nominal value. Values 
of dimensionless stiffness could be higher than the nominal value by 22 percent or lower 
by 23. 8 percent. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, February 26, 1975, 

505-04. 


REFERENCES 

1. Kirk, R. G. ; and Gunter, E. J. ; Effect of Support Flexibility and Damping on the 

Dynamic Response of a Single Mass Flexible Rotor in Elastic Bearings. 

NASA CR-2083, 1972. 

2. Lund, Jorgen W. : The Stability of an Elastic Rotor in Joiarnal Bearings with 

Flexible, Damped Supports. J. Appl. Mech. , ASME Trans., vol. 32, ser. E, 
no. 4, Dec. 1965, pp. 911-920. 

3. Rieger, Neville F. : Unbalance Response of an Elastic Rotor in Damped Flexible 

Bearings at Super Critical Speeds. J. Eng. for Power, ASME Trans. , vol. 93, 
ser. A, Apr. 1971, pp. 265-278. 

4. Hays, Donald F. : Squeeze Films for Rectangular Plates. J. Basic Eng. , ASME 

Trans., vol. 85, ser. D, no. 2, June 1963, pp. 243-246. 

5. Salbu, E. O. J. : Compressible Squeeze Films and Squeeze Bearings. 

J. Basic Eng. , ASME Trans. , vol. 86, ser. D, no. 2, June 1964, pp. 355-366. 


17 



6. Licht, L. ; and Cooley, J. W. : Dynamics of Externally-Pressurized Sliders with 

Incompressible and Compressible Films. J. Basic Eng. , ASME Trans. , vol. 86, 
ser. D, no. 2, June 1964, pp. 396-404. 

7. Tully, N. : Dampii^ in Externally Pressurized Gas Bearing Journals. The Engineer, 

vol. 222, Nov. 1966, pp. 794-797. 

8. Design of Gas Bearings. Vol. 1, Design Notes Mechanical Technology, Inc., 1966. 

9. Stiffler, A. K. : Analysis of the Stiffness and Dampii^ of an Inherently Compensated, 

Multiple Met, Circular Thrust Bearing. ASME Trans. , vol. 96, ser. F, no. 3, 
July 1974, pp. 329-336. 

10. Mullan, P. J. ; and Richardson, H. H. : Plane Vibration of the Inherently Compen- 

sated Gas Journal Bearing; Analysis and Comparison with Experiment. ASLE 
Trans., vol. 7, no. 3, July 1964, pp. 277-287. 

11. Stiffler, A. K. ; and Smith, D. M. : Dynamic Characteristics of an Inherently Com- 

pensated, Square, Gas Film Damper. ASME Trans., vol. 97, ser. F, no. 1, 

Jan. 1975, pp. 52-62. 

12. Thompson, W. T. ; Vibration Theory and Applications. Prentice Hall, Inc., p. 52. 

13. Gross, William A. : Gas Film Lubrication. John Wiley and Sons, Inc., 1962, 

p. 307. 

14. Cunningham, Robert E. ; Gunter, Edgar J. , Jr. ; and Fleming, David P. : Design 

of an Oil Squeeze Film Damper Bearing for a Multimass Flexible -Rotor Bearing 
System. NASA TN D-7892, 1975. 

15. Barrett, L. E. ; and Gunter, E. J. : Steady State and Transient Analysis of a Squeeze 

Film Damper for Rotor Stability. ME-4040-114-73U, University of Virginia, 1973. 


18 









X-Y-Y' plotter 



Figure 6. - Instrumentation schematic for air squeeze-film damping studies. 
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(a) Turbine frequency, 9500 cycles per minute; phase 
angle from meter, 90° 



(b) Turbine frequency, 28 700 cycles per minute; 
phase angle from meter, 175° 



(c) Turbine frequency, 30 500 cycles per minute; 
phase angle from meter. 177° 


Figure 7. - Phase relationship of position of exciting 
force and housing displacement. 
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HOUSING PEAK TO PEAK AMPLITUDE, IN. 




.0030|— 
.0025 — 
.0020 — 

.0015 — 
.0010 — 
.0005 — 
0 — 



TURBINE WHEEL FREQUENCY, CYCLES/MIN 


(a) Film thickness, h, 63. 5 micrometers (b) Film thickness, hv 25 micrometers 

(0.0025 in.). (0.001 in.). 


Figure 8. - Housing peak to peak amplitude and phase angle as functions of turbine wheet frequency. Supply 
pressure ratio, Pj/Pg, 2. 
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OIMENSIONIISS DAMPING. 




Ic) Restrictor coefficient. A, 14; film thickness, h, 25 micrometers (d) Restrictor coefficient, A, 6.25; film thickness, h. 38. 1 micrometers 

(0.0015 in.). * (0.001 in.). 


Figure 9. - Dimensionless damping as function of squeeze number for single pad. Supply pressure ratio, P^/Pg, 2; ratio of center region length to total damper pad length. A. 
0 . 6 . 


oi 



DIMENSIONLESS DAMPING, 


IS3 
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1 2 4 6 8 10 20 40 60 80 100 1 2 4 6 8 10 20 40 60 80 100 

SQUEEZE NUMBER, o = 12m)L2/h2pa 


(c) Restrictor coefficient, A, 4. 17; film thickness, h, 38. 1 micrometers (d) Restrictor coefficient, A, 9. 35; film thickness, h, 25 micrometers 

(0.0015 in.). (0.001 in.). 

Figure 10. - Dimensionless damping as function of squeeze number for single pad. Supply pressure ratio. P^/Pg, 3; ratio of center region length to total damper pad length, A, 

0 . 6 . 




INITIAL CLEARANCE, h, fim 

I I I ^ ^ ^ 

.8 1,2 1.6 ZO 2.4x10-3 

INITIAL CLEARANCE, h. IN. 

Figure 11. - Squeeze film damping as function of clearance. Supply pres- 
sure ratio, Ps/Pg, 2; squeeze number, o, 10. 



.8 1.2 1,6 2.0 2.4x10-3 

INITIAL CLEARANCE, h, IN. 


Figure 12. - Squeeze film damping as function of clearance. Supply pres- 
sure ratio, Pj/Pg. squeeze number, a, 10. 
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DIMENS lONlfSS STIFFNESS. K = Kh/L^Pa [ Pj/Pg) - 1 






(a) Restrictor coefficient. A, ^25; film thickness, h, 63. 5 micrometers (b) Restrictor coefficient. A, 3.51; film thickness, h, 50.8 micrometers 

(0.0025 in.). (0.002 in.). 



(c) Restrictor coefficient, A, 6.25; film thickness, h, 38. 1 micrometers (d) Restrictor coefficient, A, 14.0; film thickness, h, 25 micrometers 

(0.0015 in.). (0.001 in.). 


Figure 13. - Dimensionless stiffness as function of squeeze number for single pad. Supply pressure ratio, PJP-, 2; ratio of center region length to total damper pad length. 
A. 0.6. 



DIMENSIONLESS STIFFNESS, K = KhfL^Pg IjPg/Pa) - l' 
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(a) Restrictor coefficient, A. 1.5; film thickness, h, 63.5 micrometers 
(0.0025 in.). 


(b) Restrictor coefficient, A, 2.34; film thickness, h, 50. 8 micrometers 
(0.002 in. ). 
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(0.0015 in.). (0.001 in.). 


Figure 14. - Dimensionless stiffness as function of squeeze number for single pad. Supply pressure ratio, P^/Pg, 3; ratio of center region length to total damper pad length, A, 
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STIFFNESS COEFFICIENT, K, LB/IN. 



INITIAL Ci£ARANCE, h, \xm 



INITIAL CLEARANCE, h, IN, 

Figure 15. - Stiffness coefficient as function of initial clearance at 
squeeze number a of 10. 
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